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We calculate the next-to-leading-order (NLO) QCD correction to e + e~ — > J/ipcc at the B facto- 
ries, and present theoretical predictions on the momentum and production angular distribution for 
J/ip production, and momentum distribution for J/ip polarization at NLO for the first time. By ap- 
plying Brodsky, Lepage and Mackenzie scale setting for the renormalization scale, it is found that the 
QCD perturbative expansion is significantly improved with the unique scale choice fj* — 1.65GeV. 
Together with the ip' feed-down contribution, the total cross section and momentum distribution 
can account for the recent experimental measurement by the Belle collaboration. The total cross 
section and momentum distribution are also found to be consistent with the experimental measure- 
ment in the previous study on e + e~ — > J/ipgg. However, the production angular distribution of 
J/ip production for either the J/ipcc or the J/ipgg channel has a quite different shape in contrast 
with the new experimental data, although it fits with the experimental data when the two channels 
are added together. This situation is difficult to explain. To clarify the puzzle of J/ip polarization, 
further experimental measurements are strongly expected to testify our predictions on the momen- 
tum distribution for J/ip polarization. Our total cross section agrees with that given in the previous 
study of Zhang and Chao by using their renormalization scheme and input parameters. 
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PACS numbers: 12.38.Bx, 13. 66. Be, 14.40.Gx 



I. INTRODUCTION 

For heavy-quarkonium production and decay, a naive 
perturbative QCD and nonrelativistic factorization treat- 
ment is applied straightforwardly. It is called color- 
singlet mechanism (CSM) To describe the huge 
discrepancy of the high-p t J/ip production between the 
theoretical calculation based on CSM and the experi- 
mental measurement at Tevatron, a color-octet mecha- 
nism Q wasproposed based on the nonrelativistic QCD 
(NRQCD) 3]. In the application, J/tp related produc- 
tions or decays are very good places for two reasons, theo- 
retically charm quark is thought to be heavy enough and 
charmonium can be treated within the NRQCD frame- 
work, experimentally there is a very clear signal to de- 
tect J/tp. Now the integrated luminosity is more than 
850 fb -1 at the Belle detector at the KEKB and it is 
about 20 times larger than the integrated luminosity 
32.4 fb -1 , based on which many J/ip production pro- 
cesses was observed 0, H, II] ■ Therefore it supplies a very 
important chance to perform systematical study on J/ip 
production both theoretically and experimentally. A re- 
cent review on the situation can be found in Ref . Q • 

The measurements for exclusive J/ip productions 
e+e~ — > J/iprj c , J/ip J/ip, J/ipXcJ at the B factories 
have shown that there are large discrepancies between the 
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leading-order (LO) theoretical predictions [1, [i| 
NRQCD and the experimental measurements 
It seems that such discrepancies can be resolved by in- 
troducing higher order corrections [III [H, 03, EH, E, E3, 
liH Il9| : next-to-leading-order (NLO) QCD corrections 
and relativistic corrections. 

The cross section for inclusive J/ip production in e + e~ 
annihilation was measured by BABAR 0, [2(| as 2.54 ± 



0.21 ± 0.21 pb and Belle [4 Has 1.45 ± 0.10 ± 0.13 pb. 
Many theoretical studies [alii HI US US EE 
EE Eil have been performed on this production at LO 
in NRQCD and the results for inclusive J/ip production 
cover the range 0.6 ~ 1.7 pb depending on parameter 
choices. A further analysis by Belle [5j gives 



J/iP 



X] = 0.8718;^ ± 0.17 pb, 



(1) 



and hints a[e + e~ — > J/ip + A(non cc)] ~ 0.6 pb. It 
is suggested in Ref. [3(| that different J/ip production 
mechanisms can be tested by measuring J/ip polariza- 
tion. And a study of J/psi polarization in B — > J/ip + X 
with the BELLE detector presented in Ref. [3l| . 

For the non cc part, the contributions from the color- 
singlet channel e + e _ — * J/ipgg and color-octet chan- 
nel e + e~ — > J/ipg are about 0.2 and 0.27 pb, respec- 
tively, at the LO in NRQCD However, the signal of 
the color-octet was not found in the experiment [J, [2(| • 
Therefore, the experimental measurement by Belle is 
about 3 times larger than the theoretical prediction 
from the color-singlet at LO, and can be much more 
than 3 times by BABAR. The NLO QCD corrections 
to e + e~ — * J/ipgg has been studied by two individual 
groups recently [32, [33| , which boost the cross section to 
0.373 — 0.496 pb depending on parameter choices. Mean- 
while, a new measurement with higher statistics reported 
by Belle [34| gives 

cr[e+e- -> J/iP + X(noncc)} = 0.43±0.09±0.09 pb, (2) 

which fits well with color-singlet predictions at NLO. 

For the cc part, the experimental data by Belle 0], 
0-87±q*1 ± 0.17 pb, is about 5 times larger than the LO 
NRQCD prediction 0. However, this large discrepancy 
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FIG. 1: Typical Feynman diagrams at LO. 



was partially resolved by considering both NLO correc- 
tion and feed-down from higher excited states [HI]. It 
is also pointed out in Ref. [36| that the color transfer in 
associated heavy-quarkonium production may give im- 
portant contribution to this process. The recent experi- 
mental measurement [34[ gives 

a[e + e- -> J/ip + cc] = 0.74 ± 0.08±g g| pb, (3) 

which is even closer to the theoretical predictions. Since 
the calculation of the NLO QCD correction to this pro- 
cess is quite complicated and plays a very important role 
to explain the experimental data, it is desirable to have 
an independent calculation. A more important point is 
that there are already the momentum and production 
angular distribution for J/ip production obtained in the 
new measurement to be compared with theoretical pre- 
dictions. Furthermore, the transverse momentum distri- 
butions of J/tp polarization measured by CDF Collabo- 
ration (3?]] are still challenging our understanding of the 
heavy-quarkonium production mechanism even with the 
recent significant theoretical progresses [H, on the 



NLO QCD calculation. To understand the J/tp polariza- 
tion puzzle, there are also other J/ip related production 
processes calculated [ioj . and it is very helpful to study 
J/tp polarization in e + e~ — * J/ip + cc from both theo- 
retical and experimental parts. Therefore, in this paper, 
we present detailed study on the NLO QCD correction 
to e + e~ — > J/ip + cc by using our Feynman Diagram Cal- 
culation (FDC) package [41(, and give theoretical predic- 
tions on the momentum and production angular distri- 
bution for production, and momentum distribution for 
J/ip polarization at NLO for the first time. Our total 
cross section is in agreement with the previous result in 
Ref. [35l | when their renormalization scheme and input 
parameters are used. 

This paper is organized as follows. In Sec. II, we give 
the LO cross section for the process. The calculation of 
NLO QCD corrections is described in Sec. III. In Sec. IV, 
numerical results are presented. Further discussion on 
the renormalization scale choice is performed in Sec. V. 
The summary and discussions are given in Sec. VI. In the 
Appendix A, a trick which will bring better convergence 
in numerical calculation for some results is introduced. 
II. LO CROSS SECTION 

The related Feynman diagrams which contribute to 
the LO amplitude of the process e + (pi) + e~(p 2 ) — > 
J/^iPz) + c {pa) + c(ps) are shown in Fig. [IJ while the 
others can be obtained by reversing the arrows of quark 
lines. In the nonrelativistic limit, using the NRQCD fac- 
torization formalism, the differential cross section is ob- 
tained in n — 4 — 2e dimensions as 
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where s is the squared center-of-mass energy, e c and 
m c are the electric charge and mass of the charm quark, 
respectively. The dimcnsionlcss kinematic variables are 



defined as 

4to 2 ' ' 4to 2 

Ai = \J s C c - 1, A 2 = 

A3 = A 1 / 2 (S, S C - C , 1) EE y/(S - S C - - l) 2 - AS C - C . (5) 
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i? s (0) is the radial wave function at the origin of J/ip. 
The approximation Mj/^ = 2m c is taken. Our results at 
LO are consistent with those in Refs. [1, [35| 



III. NLO CROSS SECTION 

At NLO in a s , there are virtual corrections from loop 
diagrams. Dimensional regularization has been adopted 
for isolating the ultraviolet (UV) and infrared (IR) sin- 
gularities. UV-divergences from self-energy and trian- 
gle diagrams are canceled upon the renormalization of 
QCD. Here we adopt same renormalization scheme as 
Ref. (l3| . The renormalization constants Z m , Z2 and 
Z3, which correspond to charm quark mass m c , charm 
field tp c , and gluon field are defined in the on-mass- 
shell (OS) scheme while Z g for the QCD gauge coupling 
a s is defined in the modified-minimal-subtraction(MS) 
scheme: 
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Although the Feynman diagrams are similar, the cal- 
culation of tensor and scalar integrals is much more com- 
plicated than that in Ref. [HI]. Again, the calculation 
was done automatically with our FDC package 41] . 

The real corrections arise from a real gluon emission 
process, e + e~ — > J/ipccg. There are two types of dia- 
grams in this process, as shown in Fig. [3J Usually phase 
space integration for real correction processes will gener- 
ate IR singularities, which is either soft or collincar and 
can be conveniently isolated by slicing the phase space 
into different regions. We adopt the two-cutoff phase 
space slicing method [42| to decompose the phase space 
into three parts by introducing two small cutoffs, S s and 
S c . In this process, there is no collinear singularities and 
only the cutoff S s is needed. Then the real cross section 
is written as 



MC 



(8) 



where a s from the soft regions contains soft singularities 
and is calculated analytically under soft approximation. 
It is easy to find that the soft singularities for a gluon 
emitted from the charm quark pair in the S-wave color 
singlet J/ip are canceled by each other. We have 
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Two different frames are used to realize the division of 
the soft and hard noncollinear parts. The first one is in 
the center-of-mass (CM) frame of initial state particles 
e + e~ , where A\ and Aq are obtained as 



where /i is the renormalization scale, je is Euler's con- 
stant, (3q = ^-Ca — ^Tprif is the one- loop coefficient 
of the QCD beta function and rif is the number of ac- 
tive quark flavors. There are three massless light quarks 
u,d,s, and one heavy quark c, so n/=4. In SU(3) C , 
color factors are given by Tf — \-,Cf — §,Ca = 3. 
And $ = fa + (4/3)T F = (11/3)CU - (4/3)T F n ;/ where 
nif = rif — 1 = 3 is the number of light quarks flavors. 
Actually in the NLO total amplitude level, the terms 
proportion to SZ3 cancel each other, thus the result is 
independent of renormalization scheme of the gluon field. 

After having fixed our renormalization scheme and 
omitting diagrams that do not contribute, including 
counter-term diagrams, there are 80 NLO diagrams re- 
maining, which are shown in Fig. [5] Diagrams of groups 
(61) and (62) that have a virtual gluon line connected 
with the charm quark pair in J /tp lead to Coulomb sin- 
gularity ~ 7r 2 /w, which can be isolated by introducing a 
small relative velocity v = \p c —pc\ and mapped into the 
cc wave function: 
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where (3 c (c) is the ratio of momentum to energy of c(c) 
and 9 cS is the angel between the cc pair, kp and A are 
defined as 
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FIG. 2: Typical Feynman diagrams for virtual corrections. Groups (ai) — (04) are the counter-term diagrams, including 
corresponding loop diagrams. More diagrams can be obtained by reversing the arrows of the quark lines, and exchanging the 
places of the J /if) and open charm pairs in groups (a) and (b). 




FIG. 3: Typical Feynman diagrams for real correction process 
e + e~ -— > J/ifjccg. (a): Diagrams obtained by adding a gluon 
to LO diagrams. This type contains 24 diagrams; (b): Dia- 
grams with two quark lines. This type contains six diagrams. 



The other way is to do the calculation in the CM frame of 
the open cc pair, and it leads to much simpler expressions 
for Ax and Aq as 



Ac = 
where 
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is the ratio of momentum to energy for c or c in the CM 
frame of cc. The hard noncollinear part a HC is IR finite 
and can be numerically computed using standard Monte- 
Carlo integration techniques. The expressions for ^0,1 in 
the CM frame of cc in Eq. (fT5)) can also be found in 
Ref. 42[ . The real cross section a R is frame independent 
and should be the same no matter in the CM frame of 
e + e~ or cc. It is obviously checked in our numerical 
calculations. 



Finally, all the IR singularities are canceled analyt- 
ically. After adding all the contribution together, the 
cross section at NLO is expressed as 
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where (3q is the one-loop coefficient of the QCD beta func- 
tion. 



m c (GeV) 


a a {p) 


<7^ U »(pb) 


a(s) 


<j (1) (pb) 




1.4 


0.267 


0.224 


8.19 


0.380 


1.70 


1.5 


0.259 


0.171 


8.94 


0.298 


1.74 


1.6 


0.252 


0.129 


9.74 


0.230 


1.78 



TABLE I: Cross sections with different charm quark mass m c 
with the renormalization scale fi = 2m c and */s = 10.6 GeV. 
a(s) is defined in Eq. (|16|l . 



To study the polarization of J/ip production, we define 
the angular distribution coefficient A as the a in Eq. (2.1) 
of Ref. 

S(Pj^)[1 + A(Pj^)cob0\, (17) 



d cos 6dP 



and the polarization factor a is defined as 

da T /dPj /4 , - 2da L /dP J/i , 



dar/dPj/^ - 2da L /dPj/^ ' 



(18) 



where Pj/ip and are the 3-momentum and production 
angle of J/ip in the laboratory frame. <jt and <jl are 
the transverse and longitudinal polarized cross section. 
To calculate a(Pj/^), we use the same method to rep- 
resent the polarized cross section as Eqs. (8) and (9) in 
Ref. . This method is found numerically unstable in 
a small region of phase space due to the cancellation of 
large numbers. Therefore, the momentum distributions 
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FIG. 4: Cross sections of e + e~ — > J/ipcc as a function of the 
renormalization scale /i. The mass of charm quark is chosen 
as 1.4, 1.5 and 1.6 GeV, respectively. 

for A and a contain potentially large numerical errors in 
our calculation for Pjm < 1 GeV or Pj/^p > 4.2 GeV. 
As regards the total cross section and momentum distri- 
bution of J/tp production, a simplified method (see more 
details in Appendix. [AJ is used to calculate the amplitude 
square with very good convergence behavior in numerical 
calculations. But it cannot be applied to the calculation 
of A, a and cos# distributions. 

In the NLO calculation, we should adopt a s in the 
two-loop formula 

q.0*) = i ft ^H^/^qcd) (m 

At: /3 ln( M 2 /A 2 QCD ) /3 3 lnV/A 2 ^) ' 
with number of active quark flavors nt = 4 and A^^- = 

n J MS 

0.338 GeV. The value of the wave function at the origin 
of J/ip is extracted from the leptonic decay widths: 



IV. NUMERICAL RESULTS 

By using r ee = 5.55 KeV, together with a = 1/137, 
M J/V , = 2m c = 3.0 GeV and a s = 0.26, |i? s (0)| 2 = 
0.944 GeV 3 is obtained. For other value of m c , it should 
be multiplied by (m c /1.5 GeV) 2 . The numerical results 
are shown in Table. |TJ They are a bit smaller than those 
given by Zhang and Chao in Ref. [35[. It's just because 
of the differences in the value of R s (0) and renormal- 
ization scheme. If we choose their value of R s (0) and 
renormalization scheme, both calculations reach a very 
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FIG. 5: Cross sections for J/ipcc, J/ipgg and total, as a func- 
tion of the renormalization scale jj,. 
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FIG. 6: Cross sections for J/tpcc, J/ipgg and total, as a func- 
tion of the center-of-mass energy of e + e~ y^s with the renor- 
malization scale /i = \/s/2. 



good agreement. Thus our calculation confirms the re- 
sults in Ref. (35|. Thereafter, if not specified, we use 
a/s = 10.6 GeV, m c = 1.5 GeV and /i = 2m c as our de- 
fault choices in our results presented below. All of the 
results for J/ipgg are from our previous work [32J. 

The scale dependence of the cross section is shown in 
Figs. |4] and [5] We can see from Figs.|4]that the scale de- 
pendence of the total cross for J/ipcc has not improved at 
NLO. The curves marked with "total" in Fig. denotes 
the combination of J/ipcc and J/ipgg channels, together 
with the contribution from the feed-down of ip' by multi- 
plying a factor of 1.29. The treatment is applied to all the 
"total" results throughout this paper. Figs. [S] and [7] show 
the \fs dependence of the cross section, with /i = i/s /2 
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FIG. 7: Cross sections for J/ipcc, J/ipgg and total, as a func- 
tion of the center-of-mass energy of e + e~ */s with the renor- 
malization scale /i = 2m c . 



FIG. 8: Momentum distribution of J/ip production. The com- 
parison with experimental data will be made later. 



and /x = 2m c , respectively. We see that the cross sec- 
tion of the J/ipcc changes much milder than that of the 
J/ipgg channel as the center-of-mass energy increases. 
The asymptotic behavior of the LO total cross section in 
the threshold region for both channels can be obtained 
easily as 
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_ total(NLO) 

... total(LO) 

. . J/> + gg(NLO) 

J/V+gg(LO) 

_ J/i* + cc(NLO) 

... J/tA + cc(LO) 




with £ cE = yj — 2 and ^ gg — sfh — 1. It is clearly shown 
that the production of J/ipcc is through the p-wave chan- 
nel and that of J/tjjgg is through the s-wave channel near 
the threshold region. A simple J PC conservation anal- 
ysis can easily explain these behavior. Even more, the 
threshold of the J/ipcc channel is at y/s — Am c while 
that of the J/ipgg channel is at y/s = 2m c . 

The momentum and angular distributions of J/ip pro- 
duction are shown in Figs. [5] and [5] It is found that 
the endpoint behavior of momentum distribution for the 
J/ip + gg channel was obviously changed from LO to 
NLO due to a logarithm divergent term appearing in 
the NLO calculation. We find that the shape of mo- 
mentum distributions are similar with the recent exper- 
imental data [lij], but the angular distributions are very 
different. Neither the J/ipcc nor the J/ipgg channel can 
fit the experimental data. In Figs. [TU] and QTJ the mo- 
mentum distributions of the polarization factor a and the 
angular distribution coefficient A of J /ip are shown, and 
it can be seen that there is only a slight change for both 
a and A from LO to NLO. 
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COST? 

FIG. 9: Angular distribution of J/ip production. The com- 
parison with experimental data will be made later. 



V. MORE ON THE RENORMALIZATION 
SCALE CHOICE 

One possible way, although debatable, to choose the 
renormalization scale is by following the procedures of 
Brodsky, Lepage and Mackenzie (BLM) scale setting [43[ ; 
a unique scale choice is obtained and the cross section 
at NLO is expressed as 

a« = CT (°V)[l + ^!l&(s)]. (22) 

From the relevant results listed in Table. [Til we can see 
that the K factors become smaller and the convergence 
for QCD perturbative expansions is improved. In the 
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FIG. 10: Momentum distribution of the polarization param- 
eter a of J/ip. 
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FIG. 11: Momentum distribution of the angular distribution 
parameter A of J/i/'- 



case of m c — 1.4 GeV, the total cross section is 0.381 pb 
at LO and 0.540 pb at NLO with the K factor 1.42. To 
further consider the contribution from ip' feed-down by 
introducing a factor 1.29, the cross section is 0.70 pb at 
NLO, To include the contribution from e + e~ — * 2j* — > 
J/ipcc given in Ref |44|, we should add about 0.03 pb 
and the total cross section is 0.73 pb now. The result 
can well explain the recent experimental measurement 
0.74 ± 0.08±o:o8 pb given by the Belle collaboration [H. 
It should be mentioned that the optimal scale choice n* = 
1.65 GeV is close to Mj/^,/2, half of the typical hard scale 
in the process. 

In Figs. [T2"l and 1131 the momentum and angular distri- 
butions of inclusive J/ip are shown again, while this time 



FIG. 12: Momentum distribution of inclusive J/tp produc- 
tion with fj, — fi* and m c = 1-4 GeV is taken for the J/ipcc 
channel. The contribution from the feed-down of ip' has been 
added to all curves by multiplying a factor of 1.29. 
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FIG. 13: Angular distribution of inclusive J/ip production 
with /i = fi* and m c = 1.4 GeV is taken for the J/xpcc chan- 
nel. The contribution from the feed-down of tp' has been 
added to all curves by multiplying a factor of 1.29. 



m c (GeV) 


a 3 {n*) 


a^(pb) 


6(1) 


(7 (1, (pb) 




At* (GeV) 


1.4 


0.348 


0.381 


3.77 


0.540 


1.42 


1.65 


1.5 


0.339 


0.293 


4.31 


0.429 


1.47 


1.72 


1.6 


0.332 


0.222 


4.90 


0.337 


1.52 


1.79 



TABLE II: Cross sections with different charm quark mass 
m c . The renormalization scale fi — [i* is obtained by using 
BLM scale setting [43|, and b(s) is defined in Eq. ([22]) 
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FIG. 14: Momentum and angular distributions of inclusive 
J/ip production. The contribution from the feed-down of ip' 
has been added to all curves by multiplying a factor of 1.29. 



/i = fi* and m c = 1.4 GeV is taken for the J/ipcc chan- 
nel. It is clearly shown in Fig. [12] that the momentum 
distribution roughly fits the experimental data and for 
the J/ipgg part, the measurements at Pj/^, — 2.75 and 
3.75GeV are not consistent with the theoretical calcula- 
tions within the experimental error. For angular distri- 
butions shown in Fig. I13| we see that the predication for 
the total angular distribution agrees rather well with ex- 
perimental measurement, but neither the J/ipcc nor the 
J/tp non(cc) channel can fit the experimental measure- 
ments. In Fig. [TH a band of each line given in Figs. [T^] 
and [T3] is shown. The bands are obtained by verifying 
the renormalization scale and charm quark mass used in 
the calculation. For the J/ipcc channel, the boundaries 
of bands are given by the choices [i = fi*, m c = 1.4 GeV 
and /x = 4m c , m c = 1.5 GeV. For the J/ipgg part, they 
are determined by m c < \i < Am c with m c = 1.5 GeV. 
We see that the situation for the two channels is improved 
with these theoretical uncertainty bands when compared 
with experimental measurements. 



VI. SUMMARY AND DISCUSSIONS 

We calculate the next-to-leading-order (NLO) QCD 
correction to e + e~ — > J/ipcc at the B factories, and 
present theoretical predictions on the momentum and 
production angular distribution for J/ip production, and 
momentum distribution for J / ip polarization at the NLO 
for the first time. It increases the cross section from 0.171 
pb to 0.298 pb with a K factor of about 1.74 for the de- 
fault choice m c = 1.5 GeV and fi = 2m c . By considering 
its dependence on the charm quark mass and renormal- 
ization scale with /i = 2m c , the NLO cross section ranges 
from 0.230 to 0.380 pb. Furthermore, it will be enhanced 
by another factor of about 1.29 when the feed-down from 
ip' is considered. The total cross section agrees with that 
given by Zhang and Chao [35j when their renormaliza- 



tion scheme and input parameters are chosen. To further 
discuss the renormalization scale dependence, we applied 
the BLM scale setting [43[ for the renormalization scale 
and find that it improves the QCD perturbative expan- 
sion with the unique scale choice fi* = 1.65 GeV and 
the K factor changes from 1.70 to 1.42 for m c = 1.4 
GeV. Together with ip' feed-down contribution, the total 
cross section (0.73pb) and momentum distribution can 
account for the recent experimental measurement [3~4| 
when m c = 1.4 GeV and fi* = 1.65 GeV are used. 
The total cross section and momentum distribution for 
the e + e~ — * J/ipgg channel are also found to be con- 
sistent with the experimental measurement in previous 
studies [H, However, the production angular dis- 

tribution for J/ip production for either the J / ipcc or the 
J/ipgg channel has a quite different shape in contrast 
with the new experimental data, although it agrees with 
the experimental data when these two channels are added 
together. This situation is difficult to explain. 

To cut down the background from the very large elec- 
tromagnetic J/ip production 45j , more than four charged 
tracks are required in Belle's measurement. The mea- 
sured value for the J/ipgg channel is smaller than the 
real value and the correction to this effects is difficult 
to perform. This incomplete measurement could intro- 
duce uncertainty not only for J/ipgg total cross section 
but also for angular distribution and momentum distri- 
bution. Therefore, the systematic errors for J/ipgg mea- 
surements were underestimated, and it is still impossi- 
ble to conclude that the color-octet contributions from 
e+e" -> J/ip^CSo, 3 Pj) + g [13 are ruled out since 
its momentum distribution can be changed under the re- 
summation method by introducing a shape function [46| . 
To improve the measurement, it is better to cut the J/ip 
momentum at the large end point to cut down the large 
electromagnetic J/ip production, then the measurement 
will be a complete one without needing any correction. 
The better way to confirm or rule out the color-octet pre- 
diction is to perform the measurement in the way sug- 
gested in Ref. [Hj]. 

To clarify the puzzle of polarization in J / ip production 
at the hadron collider, a detailed study on J/ip polar- 
ization at e + e" will be very helpful. Therefore, further 
experimental measurements are strongly expected to tes- 
tify our predictions on the momentum distribution for 
J/ip polarization. 

We thank Y. J. Zhang for helpful discussions. This 
work was supported by the National Natural Sci- 
ence Foundation of China (No. 10775141) and Chinese 
Academy of Sciences under Project No. KJCX3-SYW- 
N2. 



APPENDIX A: A TRICK TO SIMPLIFY THE 
CALCULATION 

Considering the process e + (pi) + e~{p2) —> 7*(fc) — > 
H{P) + X, where X can be more than one particle, as 
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FIG. 15: Diagram of the process where the simplified method 
can be applied. 

show in Fig. [T5l Thus the matrix element is expressed 
as M = l^Mfj.. Then the square of matrix element is 
expressed as |M| 2 = L^M^M*. Now define 

M; v = J [dP x ]M M MZ, (Al) 

where [dPx] denotes the integration over the momen- 
tums of all final state particles except H. Then M' is 
expressed as 

(A2) 

The current conversion demands 

fc"M; = WM'^ = 0, (A3) 
thus M^j, is further expressed as 

M'^ = ci(^k u -k 2 gfJ , u ) (A4) 
+c 2 [(fc • P)g lia - P^k a ][(k ■ P)g va - P v k a \. 

Also, we have 

L^ v = 4:(-^k 2 g^ - -q^qu + k^q v + q^k u ), (A5) 
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